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Abstract 

A search for the Higgsstrahlung process e + e~ — * HZ is carried out, covering decays of 
the Higgs boson into any quark pair, a gluon pair or a tau pair. The analysis is based on 
the 630 pb -1 of data collected by the ALEPH detector at LEP at centre-of-mass energies 
from 189 to 209 GeV. A 95% C.L. lower mass limit of 109.1 GeV/c 2 is obtained for a 
Higgs boson cross section equal to that expected from the Standard Model if the Higgs 
boson decays exclusively into hadrons and/or taus, irrespective of the relative branching 
fractions. 
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1 Introduction 



Unlike at LEP 1 energies [Q , ALEPH searches for the Standard Model Higgs boson at LEP 2 || 
were performed under the assumption that the Higgs boson decays predominantly into bb. 
Invisible final states, which would arise for instance from a decay into a neutralino pair, were 
also investigated |2|. It is also possible, in the MSSM || as well as in more general two-Higgs- 
doublet models, to find parameter sets for which the decay into bb is strongly suppressed, to 
the benefit of other decay modes such as cc, gg or r + r~. It is shown in this letter that either 
existing or slightly modified ALEPH searches for the Higgsstrahlung process e + e~ — > HZ are 
sensitive to these decays. These analyses are based on the 630 pb -1 of data collected by ALEPH 
between 1998 and 2000 at centre-of-mass energies ranging from 189 to 209 GeV (Table [l]). 



Table 1: Integrated luminosities, centre-of-mass energy ranges and mean centre-of-mass energy values for data 
collected by the ALEPH detector from 1998 to 2000. 



Year 


Luminosity (pb L ) 


Energy range (GeV) 


(V5> ( GeV ) 


2000 


11.2 


207 - 209 


208.0 




122.6 


206 - 207 


206.6 




80.0 


204 - 206 


205.2 


1999 


45.2 




201.6 




86.3 




199.5 




79.9 




195.5 




28.9 




191.6 


1998 


176.2 




188.6 



This letter is organized as follows. A brief description of the ALEPH detector is given in 
Section 2. The event selections pertaining to a flavour-independent search for the Higgs boson 
produced via the Higgsstrahlung process are examined in turn in Section 3, and the results are 
summarized in Section 4. 

2 ALEPH Detector 

A detailed description of the ALEPH detector and its performance can be found in Refs. M 
and ||. The tracking system consists of a silicon vertex detector, a cylindrical drift chamber 
and a large time projection chamber (TPC), immersed in a 1.5 T axial magnetic field provided 
by a superconducting solenoidal coil. With these detectors, a transverse momentum resolution 
of Sp t /pt = 6 x 10~ 4 p t © 5 x 10~ 3 (p t in GeV/c) is achieved. 

An electromagnetic calorimeter placed between the TPC and the superconducting 
coil identifies electrons and photons, and measures their energies with a resolution of 
5E/E = 0.18/ \[E + 0.009 (E in GeV). The iron return yoke is instrumented with 23 layers 
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of streamer tubes and serves as a hadron calorimeter and muon filter. Two additional double 
layers of streamer tubes outside the return yoke aid the identification of muons. 

An energy flow algorithm || combines the information from the tracking detectors and the 
calorimeters and provides a list of reconstructed charged and neutral particles. The achieved 
energy resolution is cr(E) = O.QVE + 0.6 (E inGeV). 

3 Event selections 

In this analysis, decays of the Higgs boson to hadrons or to tau pairs are considered. Hadronic 
Higgs boson decays are searched for in four-jet, missing-energy and leptonic (electron or muon 
pair) final states as for the Standard Model search for the e + e~ — > HZ process. The event 
selections used here for the missing energy and leptonic final states are based on those used 
in the previous searches || [Tj, however without b-tagging information. A specially designed 
flavour-independent selection is used for the four-jet channel. Higgs boson decays to a tau pair 
are searched for in the final state r + r~ qq, using the same selection described in Ref. |7j]. 

Signal efficiencies and background contributions from Standard Model processes are 
estimated with simulated event samples which include a full simulation of the ALEPH detector. 
To study the signal efficiency for H — > hadrons, events from the HZ process are generated in 
which the H decays to bb, cc or gg and the Z into a pair of quarks, neutrinos, electrons or 
muons. Events in which H decays to a pair of taus and the Z to a pair of quarks are used 
to determine the tau channel efficiency Signal events are simulated with the Monte Carlo 
generator HZHA ||, for mn from 40 to 115GeV/c 2 in steps of 5GeV/c 2 . The simulated 
background event samples are identical to those used in Ref. ||. 

3.1 Leptonic and missing energy final states 

The leptonic channel event selection, which does not include b-tagging information, is 
unchanged with respect to Ref. ||. The reconstructed Higgs boson mass, computed as the 
mass recoiling against the lepton system, is used as a discriminant variable in the confidence 
level calculation. When this analysis is applied to the data 0, [7], fLO] , 70 events are observed, 
in agreement with the 73.4 events expected from the Standard Model backgrounds. The signal 
efficiencies for H — > bb, cc or gg are found to be quite similar, at about 80% over the entire 
mass range, except when approaching the kinematic limit for HZ production where it falls to 
40%. 

As for the leptonic analysis, the Standard Model missing energy event selection applied to 
data collected in 2000 does not include b-tagging information, and is therefore used in this 
search. In this analysis, the reconstructed Higgs boson mass is used as a discriminant variable 
in the confidence level calculation. Prior to 2000, the Standard Model missing energy event 
selections relied explicitly on b tagging, and are therefore inappropriate for this analysis. A 
modified version of the three-neural-network analysis described in Ref. is applied to this 
data sample. This analysis uses the seven-variable anti-qq and three-variable anti-WW neural 
networks used in the standard analysis. When these analyses are applied to the data, 177 events 
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are selected in the sample with 181 expected from Standard Model background processes. Cuts 
on the two neural network outputs are chosen to optimize the search sensitivity as a function 
of the Higgs boson mass. The signal efficiencies for H — > bb, cc or gg are found to be quite 
similar, and are about 40% over the entire mass range, falling to 20% near the kinematic limit. 



3.2 Final state with taus 

The search for the r + r~ qq final state of Ref. || is used here for Higgs bosons decaying into 
tau pairs. The selection efficiency is around 40%. The reconstructed Higgs boson mass is used 
as a discriminant variable in the confidence level calculation. A total of 27 candidate events is 
selected in the data, in agreement with 27.2 expected from background processes. 



3.3 Four-jet final state 

The preselection of the flavour-independent four-jet selection is similar to that of the standard 



four-jet analysis |I0| . Events are required to have at least eight charged particle tracks, and 
the total energy of the charged particles must be larger than 10% of the centre-of-mass energy. 
Events from radiative returns to the Z resonance, in which a photon escapes undetected down 
the beam pipe, are rejected by requiring the momentum p z of the event along the beam axis to 
satisfy p z < 1.5(m vis — 90), where m vis is the total visible mass in the event, expressed in GeV/c 2 . 
Events are then clustered into four jets using the Durham jet-clustering algorithm ||. The 
transition from four to three jets is required to occur for 7/34 > 0.008. Events from radiative 
returns to the Z with a photon in the detector are rejected if more than 80% of the energy of 
any jet is in the form of electrons and photons. Events from semileptonic decays of W + W~ are 
rejected by requiring that the energy of the most energetic identified electron or muon is less 
than 20 GeV. To avoid overlap with the leptonic selection, events containing a pair of identified 
electrons or muons with an invariant mass greater than 40GeV/c 2 are rejected. After this 
preselection, signal efficiencies for H — > bb, cc and gg are of the order of 70%. The numbers of 
events expected from background processes and the numbers of candidate events observed in 
the data are reported in Table |[ The comparison indicates a reasonable agreement between 
data and the expectation from Standard Model processes at the preselection level. 

This analysis uses several of the kinematic variables used in the standard analysis and in 
addition, three variables based on di-jet mass information: 



• the significance of the distance to the W W hypothesis, defined as 

£ ww = min {(m, i +m (; r2mw) 2 / (7 s 2 + m ki) 2 /^d}^ 

i,j,k,l=l,4 

where a s = 4 GeV/c 2 and ad = 10 GeV/c 2 are the resolutions on the sum and the difference 
of the di-jet masses for W + W~ production, and k, I denote the four jets reconstructed 
in the event; 
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Table 2: Numbers of expected events from background processes and numbers of candidate events collected 
at centre-of-mass energies from 189 to 209 GeV, at the preselection level for the four-jet final state selection. 



v/i (GeV) 


Background process contributions 


Data 


WW 


qq 


ZZ 


Total 


188.6 


1002.1 


261.5 


63.8 


1327.5 


1242 


191.6 


165.2 


41.9 


11.8 


218.9 


221 


195.5 


459.8 


108.1 


35.7 


603.4 


614 


199.5 


492.3 


108.0 


40.0 


640.1 


624 


201.6 


238.9 


51.2 


19.7 


309.7 


261 


204-209 


1251.0 


247.9 


102.8 


1601.6 


1601 


All v/i 


3609.3 


818.6 


273.8 


4701.7 


4563 



• the probability density functions S mii (£}iz) and £> mH (£nz) for signal and background. 
Here, £hz is the significance of the distance to the HZ hypothesis. It depends on mn and 
is defined as 



c . j [{rriij + m M ) - (m z + m H ) , [(nty - m M ) - (m H - m z ) 1 
£ H z = mm <^ J - = + J - = *-), 

i,j,k,l=l,4 I erf; 

where as and <ta are the resolutions on the sum and the difference of di-jet masses for 
HZ production. Simulated event samples are used to parametrize S mil and B mii . 



These three variables are combined with the smallest jet energy E^ n , the largest jet energy 
E™ x , and the product E^ n 8ij of the smallest jet energy and the minimum angle between any 
two jets, in a six-variable neural network. In order to optimize the performance over a wide 
range of masses, separate neural networks are trained for several Higgs boson mass hypotheses, 
ranging from 40 to 115 GeV/c 2 , in steps of 5 GeV/c 2 , at three different centre-of-mass energies : 
189, 199.5 and 206.7 GeV. For each of these networks, an optimization using the N g5 [fTl| 
prescription is performed to determine the appropriate cut for the neural network output. As 
the di-jet mass information is included in the neural network, it is not included again as a 
discriminant variable when computing the confidence level. 

As neural networks are trained every 5GeV/c 2 , a sliding method is used to determine 
the selection at intermediate Higgs boson masses. For a mass mn intermediate between two 
training masses mi and m 2 , S niH (B mH ) is interpolated from S mi (B mi ) and S m2 (B m2 ). These 
quantities are input to the network trained at ni\ and to the network trained at m^. The two 
neural networks outputs are then interpolated to calculate iViV mfl . The cut value on NN mfi , the 
signal efficiency and background expectation at the intermediate mass are similarly obtained 
by interpolations. The validity of this procedure has been established with test samples of 
signal events simulated at y/s = 199.5 GeV for masses between 60 and 100GeV/c 2 in steps 
ofl GeV/c 2 . 
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Figure 1: Expected background (solid curve) compared to collected data (dots) by the four-jet selection as 
a function of the Higgs boson mass hypothesis at centre-of-mass energies from 189 to 209 GeV. The dotted 
curve represents the contribution from the qq background and the dashed curve represents the additionnal 
contribution from the ZZ background. 



The number of events selected by this analysis in data is compared to the background 
expectation in Fig. [I] as a function of the Higgs boson mass hypothesis. The points are 
statistically correlated, as mass resolution effects usually allow the events to contribute to 
several adjacent mass bins. As a result, a deficit in one mass region can be propagated to a 
large range of mass hypotheses, as observed in the 60 and 90GeV/c 2 regions, correlated to 
the W + W~ and ZZ deficits already described in Ref. (Hi and Ref. [O . 



4 Results 

In every channel under investigation, no departure from Standard Model expectations consistent 
with the presence of a Higgs signal is observed in the data. Lower limits on the lightest scalar 
Higgs boson mass are derived as a function of £^ ad or the product of the branching fraction to 
either hadronic jets or tau pairs and of the ratio of the production cross section to the Standard 
Model production cross section. In order to obtain a flavour-independent limit for the decay 
H — > hadrons, the smallest of the signal efficiencies for H — > bb, cc and gg is used at each 
Higgs boson mass hypothesis. For a Higgsstrahlung cross section equal to that of the Standard 
Model and for 100% branching fraction to hadrons, Higgs boson masses below 110.6 GeV/c 2 
are excluded at 95% C.L., where a limit of 110.5 GeV/c 2 is expected in the absence of signal. 

When the parameter £^ ad is allowed to vary, the result of the flavour-independent search is 
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expressed as an excluded domain in the (mn, Chad) plane, as shown in Fig. 2a. Results from 
Ref. |L3| are used to exclude Higgs boson mass hypotheses below 40 GeV/c 2 . 

A similar procedure is followed to obtain a limit on the decay H — > t + t~, and the exclusion 
in the (ran, £r) plane is shown in Fig. 2b. For £ 2 = 1, a lower limit on the Higgs boson mass 
of 112.4 GeV/c 2 at 95% CL. is obtained, where a limit of 113.9 GeV/c 2 is expected in the 
absence of signal. Under the assumption that £ 2 ad + £ 2 = 1, a 109.1 GeV/c 2 lower limit on mn 
is obtained irrespective of £ 2 . 

The dominant systematic error sources, evaluated as described in Ref. 0, are included in 
the obtained limits. The finite size of the simulated event samples, the jet energy and angular 
resolutions and the uncertainties in the signal and background cross section estimations affect 
all the topologies under investigation. In the leptonic channel, lepton identification and isolation 
are additional sources of uncertainty. For the four-jet channel, systematic uncertainties due to 
differences between data and simulation in the event selection variables are taken into account 
with an event reweighting method fllJJ. The global effect of these uncertainties is to decrease 



the hadronic limit by 190MeV/c 2 , and the tau limit by 10MeV/c 2 . 



5 Conclusions 



In order to explore nonstandard Higgs scenarios, searches for Higgs bosons produced via 
Higgsstrahlung decaying to hadrons and to tau leptons were performed. The selections are 
similar to those used in previous searches, except for the search in the four-jet final state, 
where a new analysis was designed in order to cope with hadronic Higgs boson decays without a 
flavour tag. No evidence of Higgs boson production is observed in the search for either hadronic 
or tau decays in the data collected at energies between 189 and 209 GeV. For a Standard 
Model Higgsstrahlung cross section and a 100% branching fraction to hadrons, masses below 
110.6 GeV/c 2 are excluded at 95% C.L. independent of the flavour of the Higgs boson decay. 
Results on flavour-independent Higgs boson production have also been reported by the OPAL 
collaboration [|15j with lower energy data. For a Standard Model Higgsstrahlung cross section 
and a 100% branching fraction to r + r~, masses below 112.4 GeV/c 2 are excluded at 95% C.L.. 
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